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ABSTRACT 

O ^, Context. An exciting recent finding regarding scaling relations among globular clusters is the so-called 'blue tilt': clusters of the blue sub- 
I ' population follow a trend of redder colour with increasing luminosity. 

O ' Aims. In this paper we evaluate to which extent field star capture over a Hubble time changes the photometric properties of star clusters. Given 
that field stars in early type giant galaxies are very metal-rich, their capture will make blue GCs redder and may in principle explain the 'blue 
tilt'. 

Methods. We perform collisional N-body simulations to quantify the amount of field star capture occuring over a Hubble time to star clusters 
with 10^ to 10* stars. In the simulations we follow the orbits of field stars passing through a star cluster and calculate the energy change that 
the field stars experience due to gravitational interaction with cluster stars during one passage through the cluster. The capture condition is that 
their total energy after the passage is smaller than the gravitational potential at the cluster's tidal radius. By folding this with the fly-by rates 
of field stars with an assumed space density as in the solar neighbourhood and a range of velocity dispersions cr (15 to 485 kmi"'), we derive 
estimates on the mass fraction of captured field stars as a function of environment. 

Results. We find that integrated over a Hubble time, the ratio between captured field stars and total number of clusters stars is very low (< lO""*), 
• ■ even for the smallest field star velocity dispersion cr =15 km.?"'. This holds for star clusters in the mass range of both open clusters and globular 
OS ■ clusters. We furthermore show that tidal friction has a negligible effect on the energy distribution of field stars after interaction with the cluster. 
We note that field star capture at the time of cluster formation, when the cluster potential increases with time, is more efficient. However, it 
cannot explain the trend that more massive star clusters are redder. 

Conclusions. Field star capture is not a probable mechanism for creating the colour-magnitude trend of metal-poor globular clusters. 
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1. Introduction ward explaining the trend. Numerous authors have discussed 

the possibility of star cluster self enrichment, but there is a 

An exciting recent finding regarding scaling relations among ^-^^ ^^„gg conclusions as to whether and to which extent 

globular clusters is that the colours of individual globular clus- pQ^gj^ig q^s (e.g., Frank & Gisler UMl Smith \mE 

ters (GCs) of the blue sub-population are correlated with their q^^^-^ Parmentier & Gilmore |200l] Dopita & 

luminosities. This correlation is such that brighter globulars SmithUlMl Morgan & Lake[l989l Thoul et al. |2002l Recchi 

are redder (Hams et al. |2006l Mieske et al. |2006j Strader et & Danziger|2005l Prantzos & Charbonnel|2006ll. 

al. 120061 Spitler et al. |2006l Cantiello et al.mi. The ampli- presence of contaminators like stripped nuclei can 

tude of this 'blue tilt' is about 0.03 to 0.04 mag in colour per probably not give a satisfactory explanation for the trend. Given 

mag in luminosity. Assuming coeval GCs, the trend implies a Gaussian shape of the blue colour peak (Bekki et al. |2007l 

relation between mass of the GCs and the luminosity weighted p^^g ^^^^ observed magnitude range, one 

mean metallicity of its member stars. Various mechanisms have ^^^j^ ^^e "contaminators" to actually dominate the GC 

been discussed that may offer ways towards explain the trend, ^^^pj^ ^j^j^ unlikely, given the much smaller number of 

like self-enrichment (Strader et al. |2006|) or "sample contami- ^^^^^^^^ ^^^j^j expected over a Hubble time (Bekki et al. l2003l 

nation" by stripped nuclei of dwarf galaxies (Harris et al. 120061 Mieske et al | 2006l l 

Bekki et al. JOOTi). Mieske et al. ( |2006l M06 in the following) we indicate 

Self-e nrichment, e specia lly if pressure-induced (Mieske et that also the capture of field stars in giant elliptical galaxies can 

al. 120061 Parmentier |2004|), may off^er a plausible way to- !„ principle cause such colour-magnitude trends, including the 

dependence of the trend on field star density that is detected in 

Send offprint requests to: S. Mieske M06. This is because the field star population is generally much 
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redder than the blue globular clusters. Assuming a non-linear 
dependence of the capture efficiency on globular cluster mass, 
a colour-mass trend will occur. The question is whether clus- 
ters can obtain a sufficiently large population of captured field 
stars (several percent) to create a notable trend. More generally, 
it is also of interest to which extent field star capture may ex- 
plain the multiple populations detected in massive MiUcy Way 
globular clusters (e.g. Bedin et al. 120041 D'Antona et al. 120051 
Piottoetal.|2007ll. 

In this paper, we present collisional N-body simulations 
to quantify the amount of field star capture occuring over a 
Hubble time to star clusters with a time-invariant gravitational 
potential. The setup of the simulations is presented in Sect. |2] 
the results are shown in Sect. |3] The discussion in Sect. H] in- 
cludes a comparison between our results and estimates on field 
star capture in bulge globular clusters from Bica et al. (jl997D , 
which are markedly different. 



2. Setup of simulations 

We simulate the passage of field stars through a star cluster 
by means of direct A^-body simulations, using a fourth order 
Hermite scheme with individual time-steps to follow the orbits 
of field and cluster stars. For the sake of simplicity we do not 
include the tidal field of the host galaxy in the calculations. 
This is justified since interactions are restricted to the central, 
high density parts of the cluster, where the host galaxy poten- 
tial does not play a role. The host galaxy properties are folded 
in later, when capture probabilities are derived as a function of 
the ratio between tidal radius r„v/ and cluster half-mass radius 
r/,. We assume a time-invariant gravitational potential for the 
star clusters, since the typical crossing times of a few Myr are 
very small compared to mass loss time scales (e.g. Baumgardt 
& Makino|2003l Lamers et al. l2006l ). In our calculations, the 
interactions between field stars and cluster stars are calculated 
at each time step, while the cluster stars do not interact among 
themselves, but rather feel a smooth cluster potential. Allowing 
for direct interactions between the cluster stars would have in- 
creased the required computation time to prohibitively large 
values. We have therefore not included this case in the present 
study. It is in any case unlikely that direct interactions between 
cluster stars influence our results; the interaction of the field 
stars with the cluster stars happens on a crossing time, while 
interactions between cluster stars lead to orbital changes only 
on a relaxation time, which is at least a factor of 100 longer for 
the considered clusters. 

We simulate three star clusters with particle numbers of 
N=10^, 10"*, and 10^, assuming a Plummer profile for the stellar 
density distribution. Three different values of the initial radial 
velocity at infinity of the field stars Vini are simulated, namely 
Vini=0.1, 0.33, and 1 (Tduster, where (Tduster is the velocity dis- 
persion of the star cluster For each particle number and field 
star velocity, we simulate the passage of 10"* field stars. 

The distance of closest approach p{b) as a function of the 
initial impact parameter b at infinity is given by 



where is the mass of the cluster and G the gravitational con- 
stant. For the direct numerical calculation of a given field star's 
trajectory, we restrict the distance of closest passage p{b) to 
smaller than two half-mass radii of the star cluster since larger 
impact parameters do not result in any notable interaction be- 
tween cluster and field stars and can hence be excluded from 
our runs (see Fig. [U. The initial impact parameter b of each star 
is chosen randomly within the surface area that corresponds to 
p < 2 * r^. The direct numerical calculation for a given field 
star is started at a distance of 2 rtid. The velocity and position at 
this distance are calculated by analytically integrating the or- 
bit, assuming a Keplerian ellipse with parameters based on b 
and Vini. After following each individual passage, a star is con- 
sidered as captured if its total energy E after the interaction is 
smaller than the potential energy at the tidal radius: E < ^^^Ml 
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5 10 

Minimum distance/Half-mass radii 

Fig. 1. The fraction of captured field stars vs. the minimum dis- 
tance to the cluster centre for different initial stellar velocities 
and capture criteria. For the dotted and long-dashed curves, a 
tidal radius of r,,^/ - 20r;, was assumed, while the solid and 
short-dashed curves assume no tidal field of the galaxy. For star 
clusters in tidal fields, significant capture of field stars occurs 
only for impact parameters p smaller than 2r/,. 



3. Results 

3.1. Capture probabilities 

Figs. |2] and [3] show the energy distribution of field stars after 
interaction with the cluster stars. The results for clusters with 
N^IO^ to 10^ stars are directly taken from the corresponding 
N-body simulations. In Fig.[3]we also plot results for - 10^, 
where we have re-scaled the results from lower (see below). 
This was done since calculations for large numbers of cluster 
stars require exceedingly large calculation times. 
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In the units of the calculations shown in Figs. |2] and |3] the 
mean kinetic energy of the cluster stars is 0.14582, which is 
equal to the total initial energy of field stars when v„„ = Cciuster- 

Due to the encounters between cluster and field stars, the 
energy distribution of the field stars gets broadened. If relax- 
ation is responsible for the broadening, one would expect that 
the dispersion in energy of the field stars after the passage 
should be proportional to 

E ~ V fre. 

where fp^ss is the crossing time of the field stars through the 
cluster and is the clusters relaxation time. Since frei ~ 
the energy distribution should become narrower for a larger 
number of cluster stars, which is indeed observed in the runs. 
For Vini^O.lcTciustei, the width for A^=10'* is a factor of ^/9A 
smaller than for A^=10'', and a factor of V9T larger than for 
A^=10^. The corresponding factors for a scaling with the relax- 
ation time are VtTS and V8.0, which is close to the observed 
scaling. The reason for the remaining difference could be close 
enocunters between cluster and field stars, which lead to large 
energy changes and are not correctly described by relaxation. 
Since close encounters are more important for low-mass clus- 
ters, the difference with the observed scaling should become 
smaller for larger A^, which is indeed the case. 

On the other hand, the width depends only very marginally 
on the initial velocity: the width for Vini=0.1crciustei and A^^IO"* 
is only 17% broader than for Vini^l.OcTciuster for A^=10'^. Energy 
conservation requires the velocity of a field star near the half- 
mass radius to be given by v, = ■\j'^'^liustei: ^'L- ^^'^^^ the 
crossing time is inversely proportional to v», we would expect 
that the width is 22% higher for Vini=0.1cr(,iystei, which is close 
to the above value. Those well defined dependencies allow us 
to make a robust extrapolation of the capture probabilities from 
A^=10^ to A^=10^ to estimate capture rates of massive globular 
clusters. 

To obtain estimates for A^=10* we re-scaled the energy dis- 
tributions for 10^ by reducing its width by a factor of V8.3. 
This is the reduction factor expected from fiei ~ From 
the comparisons above we would expect this reduction factor 
to be a lower limit. Thus, the number of field stars below a 
certain energy that are derived from this re-scaling will be a 
slight over-estimate of the true number The mean of the en- 
ergy distributions for the various Vini was adopted identical to 
the values for A^=10^. This is possible because the mean total 
energy of stars after interaction is identical to the initial en- 
ergy. The implicit assumption that tidal friction is negligible is 
discussed in Sect. |4] In order to achieve a statistical precision 
of the extrapolation to better than 10"'^ (i.e. the inverse of the 
number of simulated stars), we fitted f-functions to the energy 
distributions of the A^=10^ cluster and re-scaled this continuous 
function to estimate the capture probability for a A^=10^ clus- 
ter. For the fitting, we took special care to not underestimate the 
wings of the energy distribution. See Fig. H] for a comparison 
between the continuous fit and the discrete energy distribution 
of field stars for the A^^IO'' cluster and Vini - O.lcrduster- 



Table 1. Field star capture probabilities as function of field star 
initial velocity and number of cluster stars A^. Probabilities refer 
to field stars with impact parameter p < 2r/,. For A^=10^ to 
10^, the numbers are directly taken from N-body simulations 
of 10000 field star passages. For A^=10^, the probabilities are 
taken from a fit to the energy distribution of A^=10^ which was 
re-scaled to the expected width for A^= 10^ (see text and Fig.|4|i. 
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In each plot of Figs.|2]and[3]we indicate the potential energy 
at rtid =5, 10, and 20xr/,. These values are typical tidal radii for 
open or globular clusters that have a half-mass radius of a few 
pc and move on orbits which are a few kpc away from the cen- 
ter of a Milky Way like galaxy. The upper part of each plot 
shows the cumulative energy distribution of the field stars after 
the interaction. The fraction of captured field stars for a given 
rtid is determined by the intersection of the cumulative distribu- 
tion and the corresponding vertical dashed line. The numerical 
values for rtid=5, 10 and 20 xr/, are given in Table [1] As can 
be seen, the capture probability is extremely low especially for 
large initial velocities ~ crdustei . However, even for small initial 
relative velocities and large assumed tidal radii the probabili- 
ties are negligible for most cases. Only for the very low mass 
case (A^=10^) of a small open cluster, the capture probabilites 
are a few percent for field stars with initial velocities smaller 

than (Truster ■ 

3.2. Capture rates 

To transform the capture probabilities into actual capture rates, 
we need to know the number of 'fly-by' field stars. These are 
those stars that over 10 Gyr approach the cluster within 2*rh 
and within the various Vini ranges. It is important to re-iterate 
that the initial velocities Vini in Table [T] are expressed in units 
of star cluster velocity dispersion crdustei - Depending on which 
mass is assumed for a single star, these ranges correspond to 
diff'erent ranges in km/s. In the following we adopt as mass of a 
single star 0.5 solar masses, which is the typical average mass 
of stars in clusters of the investigated range. The mass range of 
the investigated cases hence ranges from 0.5 * lO-' to 0.5 * 10^ 
Mq. This covers the regime from open clusters up to globular 
clusters one magnitude more massive than the mass-function 
turn-over (2*10^ M©, e.g. Jordan et al. \2QQ7i . 

For calculating the fly-by rates, we assume a typical 
globular cluster half-mass radius of r/, = 3 pc (Jordan et 
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al. 12005 1 and |2007| l . Furthermore, we adopt a field star density 
of 0.1 Lqpc^^, which is comparable to typical values in giant 
elliptical galaxies within r^ff (Romanowsky et al. |2001| l, and 
also similar to the value in the solar neighbourhood. Assuming 
a typical M/L ratio of 2.5, this corresponds to a mass density of 
0.25 Mqpc-^. 

The resulting fly-by rates are shown in Table |2] The calcu- 
lations are done for masses between 0.5* 10^ and 0.5*10^ M©, 
and for a range of relative velocity dispersions between field 
stars and cluster. The highest velocity dispersion (485 km/s) 
represents the case of a GC orbiting a giant elliptical galaxy 
like M87 (see Mieske et al. 120061 . The lowest velocity disper- 
sion (15 km/s) represents co-rotation of an open or globular 
cluster in a dynamically cold disk. The number of fly-by stars 
is given as fraction of the number of stars in the cluster For 
high field star velocity dispersions and low cluster masses, not 
a single field star with v < cr(,iys,g,- passes the cluster over a 
Hubble time. 

The capture probabilities (assuming 0.5 solar mass stars) 
convolved with the number of "fly-by s" within 10 Gyrs yields 
the mass fraction of captured field stars. These numbers are 
shown in Table|3]for the case of rtid =20xr/,. For most cases, the 
number of captured field stars is lower than 1 single star over a 
Hubble time. Only for the more massive examples 0.5*10^ to 
0.5* lO*" Mq and the cold disk case, up to a few dozen stars will 
be captured. The mass fraction of captured stars is < 2 * 10 
for all considered cases. 

This shows that field star capture over a Hubble-time will 
not change the integrated photometric parameters of a star clus- 
ter. 



4. Discussion and conclusions 



In Bica et al. (|1997| l, the capture of field stars by a globular 
cluster orbiting the MiUcy Way bulge was calculated both ana- 
lytically and by means of simulations. Those authors find that 
over 1 Gyr, the number of captured field stars by a star cluster 
of 0.5* 10^ Mq is of the order of a few to 10% that of the cluster 
stars. This relatively high fraction of captured stars is in harsh 
contrast to our findings. Why is this? 

The reason lies in a fundamentally different approach to- 
wards estimating the capture rate. In the simulations by Bica et 
al., the star cluster is placed ad hoc into a cloud of randomly 
moving field stars. Field stars which then happen to be located 
within the radius of the cluster and with relative velocity below 
escape velocity will be captured. We note that also in M06 we 
adopted this approach to analytically estimate the number of 
captured field stars. 

This approach neglects the fact that field stars feel the grav- 
itational potential of the star cluster already before they "enter" 
the cluster By definition, even a field star with a zero initial ve- 
locity at infinity will get attracted by the star cluster such that 
its kinetic energy upon reaching the cluster center is exactly 
equal to the necessary escape energy. This means that in ab- 
sence of energy exchange with star cluster stars and in absence 
of a time variance of the cluster potential, no field star will be 
captured. It is the two-body encounters that are required for 
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Fig. 4. Cumulative distribution of total energy E,o, of field stars 
after interaction with the cluster stars for Mc = 0.5 * 10^M» 
and v^O.lcTciustei. Solid (black) histogram shows the result of 
N-body simulations. The dotted line is the cumulative version 
of a f-distribution chosen to fit the simulated distribution. This 
continuous function is used to extrapolate the capture probabil- 
ities for A^=10^ clusters (see text and Table[T]i. The two vertical 
dashed lines indicate the energy required for capture for two 
different ratios between the cluster's tidal and half-mass radii. 
For the more negative energy, it is — = 10. For the less nega- 



tive energy, it is 



20. 



any field star to obtain a negative energy. The very small num- 
ber of encounters with a large energy transfer results in such 
low capture rates as presented in the present paper. 

We compare our capture probabilities with a more recent 
study by Mints, Glaschke & Spurzem ( 120071 ). Those authors 
investigate field star scattering and capture by open clusters 
that have 200, 500, and 2000 stars. They apply analytical es- 
timates as well as Monte Carlo and N-body simulations, which 
agree well with each other There are two differences between 
their approach and ours. First, Mints et al. consider head-on 
collisions, that is impact parameters p = 0. Our estimates for 
the capture probability cover the realistic range of p < 2 * r/, 
(see Fig. [T]!, and should thus yield at most equal or lower cap- 
ture probabilities. Second, Mints et al. use £ < as capture 
criterion. We impose E < , taking into account the finite 
gravitational sphere of influence of a star cluster Also this dif- 
ference decreases the capture probabilities derived by us. 

For comparison with their results, we have to consider the 
lowest number case of N=1000 in our data, and also use E <0 
as capture condition. While the Mints etal. cases of N=500 and 
N=2000 are equally distant to N=1000 in logarithmic space, 
a comparison with our data is more realistic for the N=2000 
case. This is because the capture probabilities of Mints et al. 
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Table 2. Fraction of field stars over cluster stars f(crgei^) that approach a star cluster within 10 Gyr and within 2*rh, within 
adjacent ranges around Vini indicated in column 3. = 0.1 corresponds to the range [0 : 0.18]. = 0.33 corresponds to 
[0.18 : 0.57]. '"" = 1 corresponds to [0.57 : 1]. The star cluster is assumed to have r-^ - 3pc, consist of 0.5 solar mass stars, 
and be embedded within a field star density of 0.1 Lqpc^^ . The approach fractions (or "fly-by" rates) are calculated for different 
assumed velocity dispersions of the field stars crggi^. 
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Table 3. Number ratio n of captured stars after 10 Gyrs over total number of star cluster stars for initial velocities Vini < Icrduster 
and a range of assumed field star velocity dispersions crfieid. The given values of crduster are for the case that the simulated stars 
have 0.5 solar masses. For the tidal radius we assume rtid = 20 * r/,. The table is created by multiplying the field star fly-by rates 
from Table[2]with the capture probabilities from TablefT] and adding up the figures for ratios J'™' =0.1, 0.33, and 1.0. A ratio 
of « = statistically corresponds to <0.5 captured stars. 
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will be higher than ours due to the assumption of head-on col- 
lisions. Since the capture probabilities decrease with increas- 
ing star number, the Mints et al. case for N=2000 should thus 
correspond best to our case of N=1000. For the three values 
Vini - 0.1,0.33, 1.0 cTciustei, wc find probabilities for £ < in 
our data of 0.498, 0.244, and 0.003, respectively. The respec- 
tive probabilities from the Mints et al. study are 0.48-0.5, 0.35, 
and 0.006 (the Monte Carlo simulations from Fig. 5 of their 
paper). The values agree well, lending independent support to 
the conclusions drawn in the present paper 

Baranov ( 119751 ) argue that tidal friction causes field stars to 
loose considerable amounts of energy during passage through a 
cluster, and hence become captured. Tidal friction would make 
itself note by a net energy loss of field stars after passage 
through the cluster For the simulations presented in Fig. |2] 
there is no notable effect of tidal friction. This is because the 
simulated field stars and cluster stars have the same mass. We 
have tested the influence of tidal friction by performing one 
simulation for field stars with 5 times the mass of a cluster star, 
for the case of a A^=10"' cluster and Vini = 0. Icrduster- In Fig.|5] 
the resulting energy distribution of field stars is compared with 
the same setup for equal mass field stars. The effect of tidal 



friction is notable: the energy distribution of the more massive 
field stars is skewed towards more negative values. However, 
the effect is relatively small. The amount of captured stars in- 
creases by only 20% compared to the equal mass case. We can 
therefore state that the effect of tidal friction is not important in 
our simulations. 

Investigating possible reasons for the multiple stellar popu- 
lations found in several Milky Way globular clusters, Fellhauer 
et al. (,2006) showed that during the time of cluster formation 
- i.e. when the cluster potential gets deeper with time due to 
the contraction of a gas cloud - a significant amount of field 
stars may be trapped. However, this early trapping cannot ex- 
plain the colour-magnitude trend in GCs of early-type elliptical 
galaxies. This is because the field stars captured at cluster for- 
mation would likely be more metal-poor than the newly formed 
stars in the GC. Capture at later times, when the field star pop- 
ulations have already become metal-enriched in the course of 
galaxy-mergers, is required to explain the 'blue tilt'. 

We conclude that field star capture over a Hubble-time will not 
change the integrated photometric parameters of a star cluster, 
provided that the gravitational potential of the cluster changes 
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Fig. 5. Comparison of field star energy distribution after clus- 
ter interaction for the default case of equal mass for field star 
and cluster star(solid lines), and 5 times more massive field 
stars (dotted lines). The simulations are for a A^=10^ cluster 
and Vini - 0. IcTj-iuster, cf. Upper left panel of Fig. |2l 



only slowly compared to the cluster crossing time. Field star 
capture is not a probable mechanism for creating the colour- 
magnitude trend of old metal-poor globular clusters. 
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Fig. 2. Distribution of total energy E,o, of field stars after interaction with the cluster stars. Left panels: Star clusters with A^=10^ 
stars. Right panels: Star clusters with A^=10'* stars. Bottom panels give the energy histogram, top panels give the cumulative 
stellar distribution. Indicated in each plot are the number of cluster stars and the field star initial velocities in units of the star 
cluster velocity dispersion cr(,iustei. The initial energies corresponding to the three cases of v=1.0, 0.33 and 0.1 cr are (in model 
units): 0. 146, 0.0159, and 0.00146. They coincide with the peaks of the respective energy distributions. The requirement for field 



star capture is Etor < 



-GM, 



with rtid being the tidal radius of the cluster in the gravitational field of the host galaxy. The vertical 



dashed lines indicate different assumed ratios of where r;, is the cluster's half-mass radius. 
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Fig. 3. Energy distributions of field stars after interaction like in Fig. |2l but here for clusters with A^=10^ (left) and 10^ (right) 
stars. The curves for A^=10^ result from re-scaling the distribution of A^=10^ (see text). For high particle numbers, only a very 
small fraction of incoming field stars gets bound to the clusters. 



